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APPLICATION OF THE LEADING-EDGE-SUCTION ANALOGY TO 
PREDICTION OF. LONGITUDINAL LOAD DISTRIBUTION AND 
PITCHING MOMENTS FOR SHARP-EDGED DELTA WINGS 

By Melvin H. Snyder, Jr.,* and John E. Lamar 
Langley Research Center 

. .. . SUMMARY, 

The leading-edge-suction analogy of Polhamus has been used to develop the longi- 
tudinal load distribution of the vortex lift for delta wings. This distribution is shown to 
be similar in shape to that of the potential-flow longitudinal loading for delta wings having 
aspect ratios of 2 or less. The totals of the two theoretical distributions for delta wings 
with an aspect ratio near 1 are in good agreement with the experimentally determined 
loadings over the angle-of-attack range from 0° to 30°. The corresponding predicted 
pitching moments show slightly more stability than those measured, because of loss of 
lift near the wing tips. 

INTRODUCTION 

The nonlinear behavior of lift, drag, and pitching moment with angle; of attack is of 
considerable interest because of its occurrence for the class of sharp- edged highly swept 
wings desirable for use in supersonic cruise vehicles and, possibly, in hypersonic 
vehicles. Many theoretical methods have been developed to predict the nonpotential, or 
vortex-lift, behavior (refs. 1 to 6). Most of these methods have limited .applicability 
because of inherent assumptions, such as conical flow and slender wings (refs. 2 to 4), 
or because additional information is required to implement them (r-ef. 5), . However, the 
leading- edge -suction analogy, originally postulated by Polhamus in reference 6 for the 
vortex lift on delta wings, has been shown to be a versatile technique which is useful over, 
a broad range of planforms (refs. 6 to 8) and speeds (ref. 8). This analogy proposes that 
the low pressures associated with the leading-edge vortices produce a normal force 
which has the same magnitude as the potential -flow leadihg-edge-suction force that is lost 
because of the separation at the sharp leading edge. The suction analogy has been applied 
by Polhamus to predict the drag due to lift of delta wings (ref. 9) as well as arrow and 
diamond wings (ref. 7). : ! V - v . . • .•< •.,!.=;• . , . 

^Professor, Wichita State University. 




Recently, Nangia and Hancock (ref. 10) concluded, from wind-tunnel tests at angles 
of attack at which the vortex lift was a significant portion of the total lift, that the longi- 
tudinal distribution of the vortex lift should be the same as for the potential lift on a 
slender delta wing. In order for this conclusion to be consistent with the leading-edge- 
suction analogy, the longitudinal distribution of the potential-flow leading-edge suction 
must be the same as the longitudinal distribution of the potential -flow lift. To evaluate 
their conclusion, the present study utilizes the vortex-lattice method of reference 11 to 
predict the distribution of the potential -flow leading- edge suction and then, by use of the 
Polhamus analogy, predicts the longitudinal distribution of vortex lift. The longitudinal 
distributions of the theoretical vortex lift and potential lift are compared and the total 
lift distribution is compared with experimental results. Ranges of aspect ratios and 
angles of attack were studied to determine iimits of the validity of the Nangia and Hancock 
(ref. 10) conclusion. 
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SYMBOLS 


aspect ratio, 


wing span 


_b2_ 

Sref 


iocal wing span, that is, the spanwise distance from one ieading edge to the 
other at any chordwise position 

local chord 


Co 


c L 


reference chord, -|c r 


root chord 

section suction-force coefficient, 

station suction-force coefficient, 
Lift 


J_ dS 
qc dy 

2 dS 
qbj dx 


lift coefficient, 


q s ref 


C LL - C NL cos a 
Cm 


i. i . ... . . .... . ; . ; c r Pitching moment 

pitching- moment coefficient about 1 — — - — - 

2 qcS r ef 
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C N 

C NL 

C S 

CSL 

AC p 

Kp,K v 

N 


normal -force coefficient, 


N 


q s ref 


(see eq. (7)) 


longitudinal normal-loading coefficient, C^l = 2 ^ ACp d^j^j = 


9N 

9x 


suction-force coefficient, 


(see eq. (12)) 


qS re f 

longitudinal suction-force coefficient, c s (c/c r ) 
pressure-difference coefficient, ^ 0Wer ^ — V . PP . e . r 
coefficients used in Polhamus' suction analogy (see eqs. (2) and (3)) 
normal force 


S 

s ref 

V 


static pressure 


free -stream dynamic pressure 


suction force 


reference area, 


be! 


a 


free -stream velocity 

chordwise coordinate in wing plane, origin at apex and positive aft 
(see fig. 1(b)) 

spanwise coordinate, origin at center line and positive toward right wing tip 
(see fig. 1(b)) 

angle of attack, deg 


Subscripts: 


p potential flow 

v vortex flow 
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DISCUSSION AND RESULTS 


Vortex-Lift Theory 

From flight and wind-tunnel studies, the flow around wings having sharp leading 
edges has been observed to separate at all but the very lowest angles of attack. For 
wings having leading-edge sweepback, this separated flow rolls up into basically two heli- 
cal vortices, the cores of which are inboard of the leading edges and above the wing-chord 
plane. As shown in figure 1(a), the flow reattaches inboard of the vortices. Data indicate 
that when leading-edge separation occurs, the wing lift generated exceeds that predicted 
by potential-flow solutions. Some authors have postulated that the total lift could be 
thought of as being composed of the potential-flow lift and a lift term associated with the 
vortex flow, that is, Cl = Cl p + Cl )V - 

The vortex-flow lift expression developed in reference 6 employs this postulate and 
is based on the concept that the leading-edge suction force, which exists in attached flow, 

the leading-edge flow separates. This concept led to 

K v sin 2 a'Cos a (1) 


is rotated into a normal force when 


C L,p 


Cl = K p sin a cos 2 a! + 


where 



aC N,p 

3(sin a cos a) 


( 2 ) 


3Cq 

K v = — (3) 

a(sin2o:) 

and Kp and K v depend only on planform and Mach number. Comparisons of Polhamus’ 
expression with experiment presented in references 6, 8, 12, and 13 show good agreement 
over the range of a, including the cases in which a large proportion of the total lift comes 
from the vortex flow. Extensive experimental data were available from reference 14 for 
the A = 1.147 flat delta wing which was used in the present study for comparison with the 
calculated lift curves. Figure 2 shows the excellent agreement of these data with the pre- 
dicted values from equation (1). 

Longitudinal Lift Distribution 

For slender wings, it is useful to consider the load distribution in the longitudinal 
direction. This force distribution has been termed variously: chordwise (refs. 10, 15, 
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and 16), stream (ref. 10), longitudinal (ref. 10), lengthwise (ref. 10), and cross-load 
(ref. 17) distribution. In this paper, this loading is termed longitudinal loading in contra- 
distinction to chordwise loading, which will be reserved for a given spanwise position. 

The longitudinal loadings used in this paper are determined by the procedures out- 
lined by Nangia and Hancock (ref. 10) and detailed here. Normally, N is determined 
from wing pressure data by the expression 

N=q tC c C iCpdxd ^ <4) 


where le and te refer to the leading and trailing edge, respectively. However, it can 
alternately be expressed as 


h - ,c *£ 

From the geometric relationships which exist on a delta wing, it is known that 
b £ x 


(5) 


( 6 ) 


Hence, by combining the preceding equations, the results can be written in coefficient 
form as 



(7) 


where 

The vortex-lift longitudinal loading, however, is obtained by a different procedure. 
This difference is due to the unique character of the vortex lift in that, rather than being 
the result of distributed pressure over an area, it is better regarded as a force per length 
distributed along each of the leading edges. The magnitude of this force is the same as 
that of the leading- edge suction force - the distribution of which can be determined by 
many computer programs based on lifting- surface theory. The vortex -lattice computer 
program of reference 11 was employed to provide the suction distributions used herein. 
The lattice arrangement selected for representing the delta wings in this study consisted 
of 10 chordwise positions at each of 12 spanwise stations on a half-wing, as shown in 


9N 

3x 

q- 
4 2 


( 8 ) 
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figure 1(b). This choice of vortex lattice was the result of a systematic study of the 
effects of varying the number of chordwise and spanwise stations (with the maximum 
number of vortices limited to 120). Two criteria were used, the first being comparison 
of the calculated value of K v with that cited in reference 6. The second criterion was 
that the normalized CgL curve must not have negative values and should go to zero, 
with the slope approaching infinity at the apex, as indicated by reference 18. 

It is conventional to express the leading-edge suction force in terms of a section 
coefficient c s defined by 



(9) 


This section suction-force coefficient is provided by the vortex -lattice program as a 

V 

function of The suction force may also be expressed in terms of a station suction 

coefficient c': 


r c ^T qd (^) 


( 10 ) 


Because of the delta planform, Equations (9) and (10) may be combined to 

relate the two coefficients: 


c s c r b £ = c s cb (11) 

In order to compare the vortex-lift longitudinal loading with the potential -lift longitudinal 
loading, it is desirable to use the same form as equation (7): 


C S = 2 



(12) 


and 


S = C s qS ref 


(13) 


Equations (9) to (13) lead to the equation used to obtain CgL from the computer 
program: 



( 14 ) 
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The results are presented in terms of which by use of the suction analogy is 

/C NL \ /Ctt\ 

equivalent to — — and, consequently, to — — . The longitudinal vortex-loading 

\ C N / v \ C L/ V 

coefficient as defined by equations (12) and (14) is consistent with the previously defined 
potential loading. 


equivalent to 


and, consequently, to 


. The longitudinal vortex -loading 


In reference 10, Nangia and Hancock reported the results of a series of wind-tunnel 
tests conducted on flat and cambered sharp-edged delta-wing models of A = 1. With the 

Ctt ^nl 

use of experimental pressure data , which is equivalent to , was obtained by 

C L C N 

evaluation of equation (8). The longitudinal loading was compared with a similar distri- 
bution predicted by the subsonic lifting-surface theory of Taylor (ref. 16). Since the 
normalized curve of potential -lift longitudinal loading was very similar to the normalized 
curve of total lift, Nangia and Hancock concluded that the variation of vortex lift must also 
have the same shape. 


The present report is concerned with this conclusion. To determine its validity 
limits, a systematic study of the longitudinal-loading characteristics of flat delta wings 
was performed by using the vortex-lattice program of reference 11. The normalized 
potential -lift and vortex-lift curves obtained with this program for an A = 1 delta wing 
are compared with data from Nangia and Hancock (ref. 10) in figure 3. Figure 4 shows 
the potential-lift and vortex-lift longitudinal loadings for an A = 1.147 wing compared 
with data reported in reference 14. Both figures 3 and 4 indicate general agreement in 
the forms of the vortex-lift and potential-lift curves and close agreement with the experi- 
mental data for total lift. The experimental data show some slight variation in the shape 
of the longitudinal loadings with angle of attack, whereas the theoretical loadings are 
independent of a. 


Effects of Angle of Attack and Aspect Ratio 

In addition to the changes in the forms of the normalized longitudinal -loading varia- 
tion with angle of attack, a second angle -of -attack effect is shown in figure 5. An exami- 
nation of this figure shows that the potential- and vortex-lift longitudinal loadings change 
proportions with a in such a way that a theoretical total is produced which, in general, 
agrees in shape with the experimental distribution even though the two curves have 
slightly different centroids. Note, in particular, the more rapid growth of the vortex - 
lift longitudinal loading with a than that for the potential lift - a characteristic of the 
lift development in this separated flow regime. 

The variations of the potential- and vortex-lift longitudinal loadings with aspect 
ratio are presented in figure 6. The most interesting feature noted in a comparison of the 
two sets of curves ((a) and (b) parts) is that as the aspect ratio decreases, the two kinds of 
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loadings approach one another. In fact, it is easy to show from slender-wing theory that 
the loadings are identical in the limit as A — 0. However, as the aspect ratio is 
increased above say 2, the similarity in shape of the two loadings no longer exists. Thus, 
the conclusion of reference 10 holds only over a limited range of aspect ratio. 

Figure 6 also shows that the variation of longitudinal loading with aspect ratio is 
greater for the vortex lift than for the potential lift. This occurs because the vortex-lift 
distribution is obtained from the leading-edge -suction distribution, which is really asso- 
ciated with a point on the leading edge rather than along either the span or the root chord. 
Hence, it may be graphed for a delta wing in either normalized coordinate direction with- 
out changing the shape. The reason for this behavior is that the leading-edge suction 
exists at the leading edge and is not distributed along the chord or span, even though it is 
dependent upon the negative pressure rise along a normal to the edge. With delta wings 
of small aspect ratio, normal to the leading edge is almost the same as normal to the root 
chord. Hence, graphing the vortex -lift loading along the root chord is appropriate. How- 
ever, as the aspect ratio increases to some value (A = 16 illustrated in fig. 6) for which 
the sweepback angle becomes low, then the normal direction to the leading edge is nearly 
normal to the span (rather than to the chord), and hence, graphing the loading along the 
root chord is not appropriate. In fact, for the A = 16 delta wing, the vortex-lift loading 
resembles a span load distribution more than it does the potential-lift longitudinal loading. 


Pitching-Moment Prediction 

The pitching-moment coefficient for the potential and vortex lift can be predicted 
by using the appropriate longitudinal loadings in the following equation: 


(C ) -3£N_ 

\ m /p,v - 2 2 



The results of computations made with this equation are shown in figure 7 as a function 
of a. From the figure, it is seen that the vortex lift produces a stabilizing moment and 
becomes the larger contributor at the high angles of attack. The sum of these curves is 
compared in figure 8 with the experimental data of reference 17 (A = 1.0) and in figure 9 
with the experimental data of references 13 and 14 (A = 1.147). In each case, the theory 
predicts a nose-down moment of larger magnitude than is obtained from experiment. The 
difference is in the direction expected to be caused by the actual loading centroid being 
forward of the theoretical - caused by loss of lift near the model wing tips. The theo- 
retical curve predicted in reference 13 by use of the computer program described in 
reference 19 coincides with the present prediction in figure 9. 


8 



CONCLUDING REMARKS 


The leading-edge-suction analogy of Polhamus is used to develop the vortex-lift 
longitudinal load distribution for delta wings. This type of loading has been combined 
with that obtained from potential -flow concepts and the combination has been compared 
with that determined experimentally . Good agreement has been determined for delta 
wings with aspect ratios of 1.0 and 1.147 at several different angles of attack. 

The Nangia and Hancock conclusion (in C.P. No. 1129, Brit. A.R.C.) that similar 
longitudinal load distributions result for vortex-flow and potential -flow lifts has been 
verified by the suction analogy for delta wings with aspect ratios of 2 or less. 

The theoretical longitudinal load distributions have also been used to predict the 
nonlinear pitching -moment behavior of delta wings. A comparison with experiment shows 
that the theory predicts slightly greater stability than measured. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 20, 1972. 
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(a) Leading-edge separation vortices. 




(b) Vortex-lattice representation. 

Figure 1.- Separated and attached flow -field representations for the delta wing. 
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Figure 2.- Variation of Cl as measured and predicted for an A = 1.147 delta wing. 




Figure 3.- Effect of angle of attack on the theoretical and experimental normalized longitudinal loadings 

for an A = 1 delta wing. 



Data, ref. 14 
a, deg 



Figure 4.- Effect of angle of attack on the theoretical and experimental normalized longitudinal loadings 

for an A = 1.147 delta wing. 
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